SHORT PROOFS OF THE ELLIPTIC BETA INTEGRALS 



V.P. SPIRIDONOV 

Dedicated to Richard Askey on the occasion of his seventieth birthday 



Abstract. We give elementary proofs of the univariate elliptic beta integral 
with bases |g|,|p| < 1 and its multiparameter generalizations to integrals on 
the An and C„ root systems. We prove also some new unit circle multiple 
elliptic beta integrals, which are well denned for |g| = 1, and their p —* 
degenerations. 
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1. Introduction 

Hypergeometric type functions play a major role in the theory of special func- 
tions, see, e.g., pQ. The concept of elliptic hypergeometric integrals was introduced 
by the author in |l()j . where a univariate elliptic analog of the beta type integrals 
has been discovered. First multidimensional extensions of this integral were pro- 
posed by van Diejen and the author in |3] (two types of integrals for the C n root 
system). Other generalizations were introduced in |10| (another C n integral and 
three different integrals for the A n root system) and (a new A n integral which, 
being related to the series prescribed earlier to the D n root system, can also be 
considered formally as a D n integral). These are all elliptic beta integrals found 
so far. The first complete proofs of the multiparameter C n integral of [H] and A n 
integral of ^U] were obtained recently by Rains [5]. The present paper describes 
elementary proofs of these exact integration formulas modelled along the proof of 
the Askey- Wilson integral j2j given by Wilf and Zeilberger in ■ In addition, we 
construct some new multiparameter elliptic beta integrals well defined for |g| = 1 
and prove them together with their g-hypergeometric degenerations. We would like 
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to mention also that the A n (or "D n ") integral introduced by Warnaar and the 
author in IT] has been proved there using similar elementary means. 

We take two base variables q,p G C satisfying constraints \q\, \p\ < 1. The key 
Jacobi type theta function has the form 

6(z;p) = (z;p) 00 (pz- 1 ;p) 00 , 

where (a;p)oo = Hj*lo(l — a J > 0- Its main transformation properties are 

9(z- 1 ;p) = 8(pz;p) = -z- 1 6(z;p). 

The standard elliptic gamma function 

r(z; q , P )= J] 7 ; q P (1) 

- 1 *- 1 — zq>p K 

j,k=0 

is symmetric in bases q, p and satisfies equations 

F(qz; q,p) = 9(z;p)T(z; q,p), T(pz; q,p) = 0(z; q)T(z; q,p) 

together with the reflection equation T(z; q,p)T(pq/ z; q,p) = 1. We shall drop bases 
q and p from the notation T(z; q,p) and use conventions 

r(tz ± ) = r(tz, tz- 1 ), r(z ±2 ) = r(z 2 , z~ 2 ), 

m 

r(ti,...,t m ) = IJr(t r ;g,p). 

Similarly, for theta functions we assume that 

6(t 1 ,...,t m ;p)=6(t 1 ;p)---6(t m ;p), OQz^p) = te _1 ;p). 
It is convenient to introduce the exponential parametrization: 

2-7ri^ 27ri^2. 2ni^ 

q = e "a , p — e ; r = e w i , 

g = e , p — e "3 , r = e "3 , (2) 

where £^1,2,3 are some complex numbers. For \q\, \p\, \r\ < 1, the modified elliptic 
gamma function has the form 

(1-e "2gJ+V i; + 1 )(l - e ->igi+t r k ) 
G(u;lj) = 7T- — n~ — (3) 

jk=Q (1-e "2 - e ««igrj r fc+l) 

and satisfies difference equations 

G(u + wi;«) = 0(e 2 ^;p)G( W ;a;), (4) 
G(u + w 2 ;u;) = #(e 2 ™^;r)G(u;u;), (5) 
G(u + w 3 ; w) = e-^-'^Gfu; w), (6) 

„ , . W 2 li 1i tUi LO2 1 



where 



U1LO2 oji u>2 6t^2 6oji 2 
The representation .''> 

G( M ;u;) = e -" p (")r( e " 2 "^;f,p), (7) 
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where P(u) is the following polynomial of the third degree 

P(u + ±^) = U{U2 -^ UJI \ (8) 

k—l 

shows that the modified elliptic gamma function is well defined for \q\ = 1 with 
oj\/u>2 > in distinction from T(z;q,p). It is related to modular transformations 
for the standard elliptic gamma function 0]. 

We have the symmetry G(u; u>i, u>2, ^3) = G(u; u>2, ^>i, ^3) and the reflection 
equation 

3 

G{a\u)G{b;u>) = 1, a + b = ^2uj k , 

k=l 

following from the property f Q3fc=i u k — u ) — ~P{ U )- I n the limit p,r -> 0, the 
function G(u; u)) becomes reciprocal to the double sine function 

lim — = S(u;ujx,u)2) = ,<,-, — I .s°° ■ (9) 

P.r^o G{u;u) (e 27 ' m /^q;q) 00 

We use also the notation 

G{ui, . . . , u m ) = G(ui;<jj) ■ ■ ■ G(ui;u), 

S(ui, u m ) = S(ui;u) ■ ■ ■ S(ui;u). 

For lm(uji/uj2) > (or \q\ < 1) the double sine function can be defined as the 
infinite product However, it remains a well defined meromorphic function of u 
in the domain uj\/lo2 > (so that \q\ = 1), see, e.g., [fy. Its zeros are located at the 
points u = — ojiN — and poles occupy the lattice u = + N) + ^2(1 + N). 
Asymptotically, 

lim S(u;u>) = l, (10) 

Im(^),Im(^)^ + oo 

lim e^-^^Siu-uj) = 1. (11) 

Im (^r). Im (^-)^-°° 

2. The univariate elliptic beta integral 

We take z G C and five complex parameters t m ,m = 1,...,5, assume that 
\q\, \p\ < 1 and compose the kernel function 

( . 1 \ I\m=l F(t m Z ± , At^) 

h > • • • '* s) - r(^,A 8 ±)n 1 < m <.< 5 r(^)' (12) 

where A = Yim=i ^m- This function has sequences of poles which converge to zero 
along the points 

V = {t m q j p k , ^-V +1 p fe+1 }m=l,..., 5>J -,fc e N 

and diverge to infinity along their reciprocals V^ 1 . Denote as C a contour on the 
complex plane with positive orientation which separates V and V^ 1 . For instance, 
for \t m \ < 1, \pq\ < \A\, the contour C may coincide with the unit circle T. 

Theorem 1. 

. . dz 4iri . 
p[z,ti,...,t 5 )— = — — . (13) 

z (<?;<7)oo(p;p)oo 



4 



V.P. SPIRIDONOV 



Proof. The first step consists in establishing the following (/-difference equation for 
the kernel function: 

p(z, qt 1 ,t 2 ,...,t 5 ) - p(z, ti,...,t 5 ) = g(q~ 1 z, h, . . . , £ 5 ) - g(z, t 1 ,...,t 5 ), (14) 
where 

g(z, *!,...,*») = P(z, ti ^ uUe{tltj . p) e(^A Z ; P )V (15) 

In the limits p — > and subsequent £5 — > 0, this equation is reduced to that used 
in |12j for proving the Askey- Wilson integral |2}. The simple p — ► limit reduces 
equality Ijl3|l to the g-beta integral introduced by Rahman in [7j. After division of 
equation 1141) by p(z, t\, . . . ,£5) it takes the form 

0(Az,Az-i;p) f} 2 Oihtjip) 



•/2. 



(16) 



(z 



Both sides of this equality represent elliptic functions of logz (i.e., they are 
invariant under the transformation z — > pz) with equal sets of poles and their 
residues, e.g., 

hm 6> Az ; p) l.h.s. = . 2 . I I -tt- 1 r 

with the same result for the right-hand side. Therefore, the difference of the expres- 
sions on two sides of (|16fl is an elliptic function without poles, that is a constant. 
The latter constant is equal to zero since for z = t\ validity of Qltif l is evident. 
Now we integrate equation (|14[) over the variable z £ C and obtain 

I(qt 1 ,t 2 ,...,t 5 )~I(t 1 ,...,t 5 )= Qf ^-jQ ti,..., t B )^, (17) 

where I(ti, ■ ■ ■ , £5) = J c p(z, £1, . . . , ts)dz/z and q~ 1 C denotes the contour C scaled 
by q^ 1 with respect to the z = point. Function (|15|) has converging sequences 
of poles at z = {t m q i 'p k , A^ 1 q J p k+1 } and diverging ones at z = {t^q~i~ 1 p~ k , 
Aq~i~ 1 p ^} for to = 1,...,5 and j,k G N. Taking C = T, we see that for 
|£ m | < 1 and \p\ < \A\ there are no poles in the annulus 1 < \z\ < There- 
fore, we can deform g _1 T to T in (|17|l and get zero on the right-hand side yielding 
I(qti,t2, . ..,£5) = J(tj., . . . ,£5). Assuming that \p\, \q\ < \A\, we obtain by symme- 
try that I(pti,t 2 , ■ ■ ■ , £5) — I{t\, . . . , £5). Further transformations £1 — > <7 ±1 £i and 
t\ — » p ±x £i can be performed only if they keep parameters inside the annulus of 
analyticity of the function I{t\, . . . , £5). 

Temporarily, we take real p and q, p < q, such that p n ^ q k for any n,k £ N 
and assume that the arguments of t^ 1 , m = 1,...,5, and A ±x differ pairwise. 
Now we fix C to be a contour that encircles V and two cuts c\ = [ti,t±p 2 ], c 2 = 
\(pq/A)p~ 2 ,pq/A] and excludes their reciprocals. Now we scale £1 — > £i<7 fc , = 
1,2,..., and, as soon as £ig fc enters the interval [tip, £ip 2 ], we perform the scaling 
t\ — > t\p~ x which does not take parameters outside c\ or C2. In this way we obtain 
I(q : >p~ k ti,t2, ... ,£5) = I(h, ... ,£5) for all j, k G N such that g J p _fc G [l,p]. Since 
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such a set of points is dense, we conclude that I does not depend on t\ and, by 
symmetry, on all t m . 

Alternatively, we can expand I(t\, . . . , fs) = J^°^_ I n {t\, . . . , ts)p n and find that 
Inigkii ■ ■ ■ ) is) = In(tl-> ■ ■ ■ j £5) termwise, similar to the situations described earlier 
in [3| . The coefficients I n are analytic in the parameters near the t m = points 
(it is the convergence of the p-expansion that imposes constraints on the absolute 
values of parameters from below). Therefore we can iterate t\ — > qt\ scalings until 
reaching the limiting point. As a result, the coefficients /„ and the integral / itself 
do not depend on t\ and, consequently, on all t m . 

Thus / is a constant depending only on p and q. Its value, given by the right- 
hand side of Ijl3|l . is found after forcing C to cross over the poles at z = t^, picking 
up the corresponding residues, and taking the limit t 2 —> l/t±, like it was done in 
the residue calculus of 3 . After proving the integration formula for a restricted 
region of parameters, we can continue it analytically to the parameters domain 
allowed by the contour C. □ 

Remark 1. Integral (| 1 3|1 was proven in JHj with the help of the Bailey's 21P2 sum 
and an analytical continuation in a discrete set of parameters. Here we gave an 
elementary proof which did not use any g-series or integral identities and used 
analytical continuations in a minimal way. 



3. A MULTIPARAMETER C n INTEGRAL 

In order to treat the multiparameter elliptic beta integral for the C n root system, 
which was introduced by van Diejen and the author in [H] and tagged there as a 
type I integral, we take z = {z\, ■ ■ ■ , z n ) £ C" and 2n + 3 complex parameters 
t = (t\, . . . , t2n+3) and compose the kernel function 

„/•„ + . n \ TT 1 TT IIm=l Fjtmzf) Yl m =l ^A^m) f-,o\ 

p{z,t,C n )- [I T (~±-±) 11 r /.±2 A 7 ±) ff Tit tY { ' 

Ki<j<n \ i 3 > i=l i A2 i J 1 ll<m<s<2n+3 1 \ L m L s) 

where A = Yl^^i t m . We denote 

V = {t m q J p k , A- 1 q j+1 p k+1 } m=1 ..^ 2n+3t:jtkeN 

the set of points on the complex plane along which the poles of 118fl in 2, converge to 
zero; C is a positively oriented contour separating V and V^ 1 ; dz/z — Jl"=i dzi/zi. 

Theorem 2. 

,dz 2 n nl(2ni) n , x 

P (z, t; C n )— = \ ' . 19 

Proof. The kernel function satisfies a (/-difference equation, analogous to 114|l . 
p(z, qti,t 2 , . . . , t 2n +3; C n ) - p{z, t; C n ) 

n 

= (Sii^i, ■■■,q~ 1 z l , . . . ,z n ,t) - g l (z,t)) , (20) 

i=l 

where 
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Dividing by p(z,t;C n ), we obtain 

yr 9(t 1 zt:p) 2 ^ ] 3 e(AtJ 1 ;p) MfaAjp) 



i ^Asf ; p) JL A 2 ; p) ITS 3 ^ ; p ) ** ( ** 2 ; p) 

Both sides of this equality are invariant under the transformation Z\ — > pzi and have 
equal sets of poles (the poles at z\ = Zj,zJ l ,j = 2, . . . , n, and z\ = ±p k ^ 2 , k e N 
on the right-hand side cancel each other) and their residues, e.g., 

lim OtAz^m. orr. h.s.) = ^ TT TT ■ 



Therefore, the functions on two sides of (|22|l differ only by a constant independent 
on zi which is equal to zero since for z\ — t\ equality l\'22l is evident. 
Integrating equality (|20() over the variables z € C n , we obtain 



I(qt 1 ,t 2 ,...,t2n+3)-I(t)=Y,( I -f 

i=l \"'C i - 1 x(g- 1 C)xC™-' JC" 



?<(*,«) — , (23) 
z 



where I(t) = J Cn p(z, t; C n )dz/z and q 1 C is the scaled partner of C. 

Poles of functions l|21|) in Zi converge to zero along the set Zi = {t m q : >p k , 
A^qip^ 1 } and diverge for = {t^ n 1 q~ : >~ 1 p~ k , Aq~i~ 1 p~ k ~ 1 }, where m = 1, . . . , 2n- 
3, j,k <E N. For \t m \ < 1 and \p\ < \A\ the region 1 < \zi\ < does not con- 

tain poles, so that we can take C = T, deform <7 _1 T to T in ill'MI) and obtain 
I(qti,t2, ■ ■ ■ ,^2n+3) = I(t)- Repeating almost literally the analytical continuation 
procedure used in the n = 1 case, we see that / is a function of p and q only. The 
latter is found from the residue calculus of |3] to be equal to the right-hand side of 

(E|. □ 

Remark 2. Integration formula (|19fl was conjectured in j^] with many justifying 
arguments (e.g., for p — > it is reduced to a Gustafson's integral jHj). It was 
proved in |S] after a reduction to computations of fairly complicated sequences of 
determinants on dense sets of parameters. Our proof is essentially shorter being 
based on the completely elementary tools. 

Remark 3. As shown in [JJJ, integral (|19|l implies validity of another highly non- 
trivial C n integral with six parameters (tagged as a type II integral), which gen- 
eralizes the Selberg integral and its g-analogues due to Gustafson 5 . A direct 
derivation of the latter integral was given in |S] . It would be interesting to consider 
it also from the present paper point of view. 



4. A MULTIPARAMETER A n INTEGRAL 

For analyzing the type I A n elliptic beta integral proposed in |l()j . we take 
z = (zi,. . . ,z n ) G C", define z n +i y i a the relation Yi7=i z i = 1, and introduce 
2n + 3 complex parameters t — (ii, . . . , t n+ i) and s — (si, . . . , s n+ 2)- The A n 
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kernel function has the form 



n+1 rcti ixtmzr 1 ) n;=i 2 r&%) w) 



X r<ii +1 Tiziz^zr^m |J "WP (24) 
where T = H"ii ^ m an< ^ ^ = II j=i s j ■ This function has poles at the points 

Zi = {t m q 3 p\ (TS)" V+V+ 1 }, i = l,...,n, z-l 1 = Z 1 ---Z n = {siq ] P k }, 
with m = 1, . . . , n + 1, I = 1, . . . , n + 2, j, k € N, converging to zero and at 

* = {«rW*}, i = I,---.", = {t- 1 q- j p- k ,TSq- j - 1 p- k - 1 }, 

diverging to infinity. 

Theorem 3. Suppose that \t m \,\si\ < 1 and \pq\ < \TS\. Then 

dz (n + l)!(2^)" 
p(z,t,s; A n )— = — — . (25) 

Proof. The A n -kernel function satisfies a q-difference equation, analogous to (|14|l . 
p(z, qh 1 t 2l tn+i, s; A n ) - p(z, t, s; A n ) 

n 

= ^ (9*( Z 1> ■ • • ' g^ 2 *) • • • 1 Z "> S ) - 4 ) S )) ) ( 26 ) 



i=l 



where 



11 /)/■„.„-!. W\ 11 ()(t, S ;-rA Q(TTK?,-n\ ^ ' 



p(z,t,s;A n ) .}, L ,.9{ziZ, \p) fj: 0{t x sj\p) 0(T,TSzi\p) 
Using the relation 

p{. . . ,g~ 1 z i , . . .;A n ) _ r y? 9( Sm z n+ i;p) "y^ OJ^jZ^p) 9(q~ 1 TSz i ;p) 
p(...,Zi,...; A n ) ~ JU B^SraZtip) 1 = 1 ^g-^z^jp) 0(TSz n+1 ;p) 

TT ^(g~ 1 ^^ 7l .g~ 1 ^^n+i;P)g(g~ 2 ^^n+i;P)^ /jgj 

equation l|26(l can be rewritten in the form 
11 0(T5z i;P ) 0(T;p) 11 e(iia i; p) 

% — 1 J — 1 

t x e{t x TS\p) A i A efazj^p) 

' 2—i v .ar~. ~-i . ^ 11 



2, 



'+i "±_ 2 n n + 1 fli'/^r 1 -r ) W('«- 1 ^T/r 1 



2 



n+1 j=1 llj=2 <-j^n+l 



nS%" 1 ^»1i;p)«(^^i;p) 



,. w ^.i,) - n W 9 (r S2i;P ) J ■< 29 » 
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Both sides of this equality are invariant under the transformation zi — ► pz\ and 
have equal sets of poles in z\ and corresponding residues. Indeed, it is not difficult 
to check coincidence of residues at z\ = (TS) -1 and we skip it. Comparison of the 
residues at z~ +1 = TS yields the equality 

n n Q(rpc -1 ^ "+ 1 

n Z *{ ;..r nw;p) 

i=l i=l,#i [ 1 1,P) 3=2 

ri+1 

= d(*r 1 5- 1 ;p)n%- 1 t i TS;p). (30) 

3=2 

If we multiply both sides by IIi<i<j<ri z j@( z i/ z jiP)i then they become antisym- 
metric holomorphic theta functions of 2i,...,z n -i obeying A n _i symmetry and 
vanishing at z n = < n. But any such function is proportional to the given 
factor. Therefore the left-hand side of H3U|) should be independent of Zi. Setting 
Z\ = q(TS) -1 we obtain the right-hand side expression. 

Residues of the poles at z n +\ = t^q -1 , k = 2, . . . ,n + 1, sum to the expression 
proportional to 

(31) 

which is equal to zero being the TS = q/tk subcase of l|3l7|l . 

The poles at z\ — Zj,l/Zj,j = 2, . . . ,n, on the right-hand side of l|29|) cancel 
each other. Finally, comparison of the residues of the poles at z\ = z n+i yields the 
equality 

After multiplication of both sides by [] 1<!< ^ <11 Zj6{zi/ zf,p), they become holomor- 
phic A n -2 antisymmetric functions of Z2, ■ ■ ■ , z n —i (zi is considered as a parameter) 
vanishing at Zj — Z\ , z n and therefore they are proportional to each other up to a 
factor independent on Z2, • ■ ■ , z n -i- For Z2 = qzi, equality (|32[1 is true and, so, it is 
true in general. 

From the listed properties it follows that the difference of two sides of equality 
(|29[1 does not depend on z\. Setting z\ =t%, we see that this difference is equal to 
zero and Q29J1 is true in general. 

Integrating equation 126J1 over the variables z G T n , we obtain 



,t n+ i,s) - I(t,s) = VM / 




I(qti,t 2 ,. . . ,t n+ i,s) - I(t,s) = y J( / -/ \g l (z,t,s) — , (33) 

l x(g- 1 T)xT"- i JT" J z 

where I(t,s) = J T „ p(z,t,s; A n )dz/z. 

Poles of functions Qi{z, t, s) are located at 

Zi = {t m q j p k , (TS)- 1 q' P k+1 }, i=l...,n, m = 1, . . . ,n+ 1, 

z n+i = { s m9 J P J }i m = l, ...,n + 2, 

converging to zero, and 



Z; 



= { s mQ J 1 P fe }> i = 1 • • • ) n, m = 1, . . . , n + 2, 



Z. 



{t^q-'-ip-*, t-V J P , TSg-'-V* -1 }, m = 2, . . . , n + 1, 
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diverging to infinity for j, k 6 N. Taking \t m \ < \q\, m = 2, . . . , n+1 and \p\ < \TS\, 
we see that the region 1 < \zi\ < |<z| _1 does not contain poles and we can deform 
q~ 1 T to T in (|33|l yielding zero on the right-hand side. Thus, I(qt\, ti, . ■ ■ , t n +i, s) = 
I(t, s) in the taken parameter region. 

Expansion of the integral in the infinite series in small p allows the iteration of 
the scaling t\ — > pt\ termise until reaching the limiting point t± = 0. Applying 
the analytical continuation procedure similar to that used in the C n case, we see 
that I(t, s) does not depend on t m . Integral 1251) can be rewritten in the t m <-> s m 
symmetric form 

U7=i U^Ji^itrnZ' 1 , S m Zi) dz 



Hl<i<j<n+l^( Z i Z j i 2 i Z j) 

_ (n + l)\{2m) 



z 

n+2 n+2 



{Q . Q) n (p . p) n II nS»£,Atf) II ms k ). (34) 

where S — Y\^=i s m,A — rim=i AS = pq. Therefore I(t, s) does not depend 
on s m as well. The explicit form of I = I(q,p), given by the right-hand side of (|25|l . 
is found from the residue calculus described in □ 

Remark 4. Formula l|25|l was conjectured in |1(J| with various justifications. Forp — > 
it is reduced to a Gustafson's integral and its first proof was obtained in jS] with 
the determinantal technique similar to that used for treating the multiparameter 
C n integral. Again, our proof looks simpler. 

Remark 5. As shown in jTH], a combination of 1)19(1 and l(25|l integrals implies validity 
of two type II exact integration formulas (which are conventionally considered as the 
A n root system objects), having different expressions for odd and even multiplicities 
of integrations. Their direct analysis from the present point of view would be of 
interest as well. 

5. Modified integrals with a base on the unit circle 

All the elliptic beta integrals discussed so far are defined only when both bases 
lie within the unit circle \q\, \p\ < 1. Elliptic hypergeometric integrals which are well 
defined when one of the bases lies on the unit circle, say \q\ = 1, were introduced in 
jlOj. They are constructed with the help of the modified elliptic gamma function 
(PJ, 0- A unit circle analog of the type II six parameter C n elliptic beta integral 
was constructed in P). Below we describe several new unit circle elliptic beta 
integrals of type I for the root systems A n and C n . 

Theorem 4. For the C n root system we take n variables u — (u\, . . . ,u„) G C n , 
In + 3 complex parameters g — (g l7 . . . , g2n+z); an d define the kernel function 

, n . r , = TT n^r 3 G(9m ± uQ I\n^G(A~ 9m ) 
PK ,5, n) 11 G{±2uiiA±Ut) Ul< m<s <2n+ S G( 9m + 9 s ) 

x II 7vT 1 T y (35) 

where A = Y^m=i 9 m > ^ e m °dified elliptic gamma function G(u) = G(u:u>), and 
G(c ± a ± b) = G(c + a + b,c + a — b,c — a + b,c — a — b). 
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Now we suppose that Imiioxjuj-i) > and Im((jJ§ / uji) , Im(u)s / > 0, and 
Im(g m /ui3) < 0, m = 1, . . . , 2n + 3, Im{{A — uz)/^) > 0. 

Then 

P p(u,g;C n )^...^ = K^nl (36) 

with 

«= , " (g " ! f" , (37) 

and i/ie integration along the straight line connecting points — UJ3/2 and UJ3/2. 

Proof. After substituting relation (J7J in (|35(l and some straightforward computa- 
tions, we obtain 

-, n -p-r2ri+3 -pi, ± - ~\ 

n(ii n- fl ") - p-*inP(P) TT 1 TT llm=l 1 \ l ™ z i '' r 'P> 

P{ n) x<M< n r(^;r,p-)M T{*?M*e> 

TT rrt * -f sv ( } 

where Zj = e" 27 ™"^ 3 , i m = e - 2ni9m / U3 , and A = Um=i tm- Since 

du - U3 
3 2iri Zj 

the integrals in Uj become equivalent to integrals in zj over the unit circle T. Ap- 
plying C n integration formula l|19(l with C = T (it is allowed due to the constraints 
\t m \ < 1, \pq\ < \A\ which are assumed to hold) and q,p replaced by r,p, we obtain 
that the left-hand side of (|36|) is equal to K n 2 n nl, where 

W2(r;r)oo(p;p)oo 

Using the transformation property of the Dedekind ry-function 

e -# (e~ 2 ™' T - e ~ 2 ™/A = {-i T ) l / 2 e^ (e 2mT ; e 2mT ) , 

it can be shown [5J that k takes the form presented in l|37(l . □ 

Theorem 5. For the simplest A n root system integral, we take integration vari- 
ables u = (u±, . . . ,u n ) £ C" and complex parameters g = (gi, . . . , g n+ i), h = 
(h\, . . . , h n+ 2)- Denoting u n+ \ = — X^fc=i we define the kernel function 

n+l T-rn+1 n( „ _ „, \ TTn+2 



FELi G( 5m - Ui) I] "=1 G% + «,) G(H + 9i ) 



x TT I 1 f? G(F + H- hj ) 



where F = Y^m=i 9 m an< ^ H = Sj=i hj- Supposing that 



Im—>0, Im—>0, Im—>0, 

W2 U)l LU2 

Im ^ <0 , Imh <0 , Im F + H - UJl - UJ2 >0 

CO3 UJ3 UJ3 



ELLIPTIC BETA INTEGRALS 



11 



for m = 1, . . . , n + 1, j = 1, . . . , n + 2, we have 



2 p( u ,g,h;A n )^...^ = K n (n + l)\ (40) 

with the same k as in the previous theorem. 
Proof. Substituting relation Q in lj3ll|) . we obtain 

fi r(r^;f,p)n"i 2 rfe Si ;f,p) 

where Zj = e - 27rlu ^^, t m = e - 2 "f™/" 3 , Sj = e - 2 ^^ 3 , T = rimti*™ S = 
Y[j=i s j- Similar to the C n case, the integrals in Uj become equivalent to integrals 
in Zj over the unit circle T. Applying A n integration formula l|25|) (which is allowed 
since our constraints arc equivalent to \t m \ < 1, \si\ < 1, and \pq\ < \TS\) with q,p 
replaced by r,p, we obtain that the integral in (|40|l equals to K n (n+ 1)1. □ 

Theorem 6. For the unit circle analog of the type I elliptic beta integral introduced 
in [TT], we take integration variables u = (u\, . . . ,u n ) G C" and complex parameters 
g= (gi,...,g n+3 ), h= (hi,...,h n ). Denoting u n +i = - YJk=i u k, we define the 
kernel function 

S(u g h- A ) = TT G(F - u % - uj) yr yr G(F + h m - g 3 ) 

Ul9, ' " "Kiii+i^-^^-^HIi G(h m + 9j ) 
tt n:u Q{h m + Ui ) n;+ 3 G( 9j ~ up i 
x fi n:Ui g(f + ^ m - Ui) ^ .ii n+3 g(f - ft - 9k y { } 

where F = ^2^=1 9j- Supposing that Im(uii/ui2) > and Im{w^/u}{), Im{uis/u>2) > 
0, and Im(gj/uj 3 ) < 0, Im(h m /ui 3 ) < 0, Im((F + h m — u)\ — cj 2 )/^3) > for 
j = 1, . . . , n + 3, to = 1, . . . , n, we /iGroe 

/ 2 ••■/ ^ff.^A.)—- ■•—=«"(»» + 1)1 (43) 

im'i/l i/ie same n as in the previous cases. 

Proof. Substituting relation Q in H42() . we obtain 

<f(u, 5 , ft; A n ) = e-"" p (°) A(z, i, s; A n ), 

where 

n .P^ffn^f <«> 

l<i<i<n+l 1 Z J' r 'Wj = lm=l i ^™ S J' r 'W 

x 11 rr™ rcn+ z\ 11 r./r>.-i.-i.s ^ 



n:=i r(^ri ; r,p) 1 <,. < 1 fc i n+3 r^ 1 ^ 1 ; f,p) 
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with Zj = e - 27rai ^ ; Sj = e- 21 "*/" 3 , t m = e - 2wih ^/^, D = j-j«+3 g ._ Ag shown 
by Warnaar and the author in for \sj\, \t m \ < 1 and \pf\ < \Dt m \ (which is true 
by the taken assumptions) the following A n (or u D n " ) elliptic beta integral holds 

AW,^,^™. (45) 

Applying this formula in our case, we obtain that the left-hand side of 1}43}) equals 
to K n (n + 1)1 as required. □ 

6. (j- Reductions of the modified integrals 

In this section we consider q-reductions of the modified elliptic beta integrals 
described in the previous section. They correspond to the limit Im(£J3) — > oo taken 
in such a way that p, r —> 0. In this limit the modified elliptic gamma functions 
reduce to the double sine functions and we obtain some beta type integrals over 
infinite contours (for earlier examples of non-compact q-beta integrals of such type, 
see, e.g., P] and references therein). However, this is a formal limit and the resulting 
integrals need to be proved rigorously. 

Theorem 7. For the C n root system we take the notation used in formula 1, 9.5}) 
and define the kernel function 



S(±2ui,A±u z ) 
1_, ,. • „ »="i rim=l S(9m ± Ui) 



p(u,g;u;C n ) = JJ S(±Ui ± Uj) J[ n , (j 



Ill<m<a<2n+3 ^ (dm + 9 s) 

Ul:=iS(A-g m ) 



(46) 



where A — J>2m=i 9™ an< ^ ^( u ) = S(u;u>) is the double sine function. Now we 
suppose that Im(uji/uj2) > and Re{uJi/u)2) > together with 

Re(g m /w 2 ) > 0, Re((A - cji)/w 2 ) < 1 

for m = 1, . . . , 2n + 3. Then 

p(.,,;a;;C„)^i-^ = (-2rn!« (47) 
w 2 ^2 (q; g)So 

with the integration contour L = zol^K. 

Proof. First, we consider convergence of the taken integral. For this we assume 
that the integration domain is bounded to an n-dimensional cube and consider 
integration in the variable Uk = iw 2 Xk for some fixed k when the size of the cube 
goes to infinity. Due to the Uj —* —Uj symmetry, it is sufficient to consider the 
subdomain < Xj < 00, j = 1, . . . ,n. We impose the constraint Re(wi / u> 2 ) > 0, so 
that for Xk — * +00 we have Re(uk/it+>i) — > +00. From the asymptotic behavior of 
S(u) described in ljTU)l . (fXTft . we obtain 

S(±2uk,A ± Uk) _ _2(n-l)irfa;^W2/wi-25rnx k (l+W2/wi)\ 
Um=l S (9m ±U k ) 

Using the reflection formula 

S(u)S(-u) = e -™ B ^ u) {\ - e -aTi«/«»)(i _ g-airfuM^ 
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for Xk — ► +00 we find the estimate 

Yl S{±u t ± u,) = 0{e- 2m{n - 1] £*=i 

l<i<j<n 

independently on the values of other xj. Thus, on the boundary of the integration 
domain we have 

p(u, g; w; C n ) = ( e -2-(i+- 2 M) £™ =1 ), 

i.e., the kernel function decays exponentially fast and the integral converges. 

After taking the limit p — > in l|2U|l . we obtain the difference equation for our 
kernel function (scalings of t m and Zi by </ ±:L are equivalent to the shifts of g m and 
Ui by ±wi): 

p(u,gi +uj!,g 2 , ■ ■ .,g 2n +3;u;C n ) - p(u, g;cj;C n ) 

n 

= ^2 ...,«» -wi,...,tt7i, 5) - M u ,g)) , (48) 



1=1 



where 



/< (n >g ) = M g ^;C n )^ (1 _^ )(1 _^ j) 

with Zj = e 2 ™^, t k = e 27 " 9 */"' 2 , a = e 2mA/u 2 ^ g ince is symme tric in w li2 , 
we have a similar equation for shifts by 012- 

We suppose now temporarily that quasiperiods u)x,2 are real, incommensurate, 
and u>i < u>2- We impose also the constraint He(A) < w\. Integrating equation 
(14811 over Uj along the contour L, we obtain 

I(gi + (ji,g 2 ,...,g2n+3)=I(g)= / p(u, g; w\ C n )dui ■ • ■ du n . 

This follows from the observation that under taken constraints the functions fi(u, g) 
do not have poles in the strip — lo\ < m < and we can shift the contour L — u)\ 
to L. By symmetry, we have I(gi + u>2, 52, ■ ■ ■ , 52^+3) = 1(g)- After an appropriate 
finite deformation of the integration contour we can reach the equality I(g\ + joJ\ — 
kcu2,g2, ■ ■ ■ ,52n+3) = 1(g) for all j,k G N such that ju)\ — kuj 2 S [O,^]. Due to the 
denseness of this set of points, we conclude that the function 1(g) does not depend 
on gi and, so, on all g m . As a result, its explicit form, given by the left-hand side 
of lPT7|) , can be found from the residue calculus analogous to those described in |3j . 
Finally, by analyticity we can extend permitted values of parameters to the domain 
indicated in the formulation of the theorem. □ 

Theorem 8. For the simplest A n root system integral, we take the notation for 
parameters as in formula ^39j) and define the kernel function 

, , , , S(F + H + u l )U;^S(g t + h J ) 

t\U n J=\S( 9m -u i )U^S(h j + u i )S(H + gi ) 

x J] Siui-UM-m) S(F) f[ Jp+'^l.y (50) 

l<i<j<n+l j=l ^ i' 
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Now we suppose that Im(u>i/u) 2 ) > and Re(u)i/u> 2 ) > together with 

Re(g m /uj 2 ) > 0, Re(hj/u) 2 ) > 0, Re((F + H- u{]/w 2 ) < 1- 

Then 

ft u a \ dui duN I iw , im(M)~ fK-n 
/ p{u,g,h;uj;A n ) = (-1 ri + l!- — (51) 

with the integration contour L = iu; 2 M.. 

Proof. Convergence of the integral under consideration follows in the way similar 
to the previous case. 

Needed difference equation for the kernel function is obtained from the p — » 
limit of 

p(u,gi +uj 1 ,g 2 , ■ ■ .,g n+ i,h;u};A n ) - p(u, g, h;uj; A n ) 

n 

= '^2(fi(ui,...,u l -uj 1 ,.. .,un,g,h) - fi(u,g,h)) , (52) 
»=i 

where 

fi{u,9,h) _ JT 1 - h/Zj yr 1 - ZjSj (1 - ZjT/ti){l - TSt{)ti ^ 

p(u,g,h;oj;A n ) J-| ■ . 1 - Zi/zj +J- 1 - tisj (1 - T)(l - TSzijZi 

withzj = e 2 ™"^ 2 , t fe = e 2™W^ Sj _ e 2mhj/<^ ; and ^ = e 2 ^/" 3 . By symmetry, 
we have also a similar equation for the shifts by u> 2 . 

We suppose now temporarily that quasiperiods oj\ t2 are positive (say, u)% < u> 2 ) 
and incommensurate. We impose also the constraints Re(F + H) < lo 2 and 
Re(<7m) > Wit m = 2, . . . , n+ 1. Integrating equation (|52|) over along the contour 
L, we obtain I(g x + u>\, g 2 , ■ ■ ■ , g 2n +3, h) = I(g, h) = J h p(u, g, h; w; A n )dui ■ ■ ■ du n 
(the functions fi(u,g, h) do not have poles in the strip — Wi < Ui < 0). 

Similar to the C n case, for an appropriate choice of the integration contours 
we have I{g\ + ju)\ — kuj 2l g 2 , ■ ■ ■ , <?2n+3, h) = I(g,h) for all j, k £ N such that 
ju>i — kuj 2 € [0,0)2]- As a result, the function I(g,h) does not depend on g\ and, 
so, on all g m . Applying the residue calculus, analogous to that described in |10| . 
we find that I(g, h) does not depend on h as well and equals to the expression 
on the right-hand side of (|51|l . Again, analytical continuation in parameters from 
restricted values to the admissible domain completes the proof of the integral in 
question. □ 

Remark 6. The Gustafson's plain hypergeometric A n integral [5] is relevant for 
the ordinary quantum Toda chain model (an observation due to S. Kharchcv). 
Therefore, it is natural to expect that a simplified version of the integral Ij51|l 
is related to normalization conditions of the Hamiltonian eigenfunctions for the 
quantum multiparticle g-Toda chain of 0. 

We shall not consider here unit circle analogs and relevant g-degenerations for 
two other types of the A n integrals introduced in [TO] , 
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Theorem 9. For derivation of the q-reduction of integral 14$ > we take the notation 
for parameters as in and define the kernel function 




Re(g j /w 2 ) > 0, Re{h m /u 2 ) > 0, Re((F + h m - wi)/w 2 ) < 1- 

Then 



with the integration contour L = za^R- 

Convergence of the integral (|55|l is deduced analogously to the previous cases. 
Further steps of the proof of this theorem are similar to the described A n and C n 
situations. Using the equation for kernel function appearing after taking the p — > 
limit in the equation constructed in it can be shown that the left-hand side of 
(|55|l is not changed after the shifts g m — > g m + Wi,2- Such shifts can be iterated 
after an appropriate deformation of the contour of integration and this leads to 
the independence of the integral on parameters g m . Application of the appropriate 
residue calculus (the p — > partner of that used in JT]) leads to equality (|55|l . We 
skip the details for they are simple enough and repetitive to those given above for 
other root systems. 



This paper was written under the influence of Richard Askey's thread to search 
for elementary treatment of special functions. The author is indebted to J. F. van 
Diejen, S. M. Kharchev, E. M. Rains, and S. O. Warnaar for stimulating discussions. 
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